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ABSTRACT 

We present results of the analysis of near infrared spectroscopic observations of 6 high-redshift 
quasars (z > 4), emphasizing the measurement of the ultraviolet Fen/ Mgn emission line strength 
in order to estimate the beginning of intense star formation in the early universe. To investigate 
the evolution of the Fell/Mgll ratio over a wider range in cosmic time, we measured this ratio 
for composite quasar spectra which cover a redshift range of < z < 5 with nearly constant 
luminosity, as well as for those which span ~ 6 orders of magnitude in luminosity. A detailed 
comparison of the high-redshift quasar spectra with those of low-redshift quasars with comparable 
luminosity shows essentially the same Fell/Mgll emission ratios and very similar continuum and 
line spectral properties, i.e. a lack of evolution of the relative iron to magnesium abundance of the 
gas in bright quasars since z ~ 5. Current nucleosynthesis and stellar evolution models predict 
that a-elements like magnesium are produced in massive stars ending in type II SNe, while iron 
is formed predominantly in SNe of type la with intermediate mass progenitors. This results in 
an iron enrichment delay of ~ 0.2to0.6Gyr. We conclude that intense star formation activity 
in the host galaxies of z > 4 quasars must have started already at an epoch corresponding to 
Zf ~ 6 — 9, when the age of the universe was ~ 0.5 Gyr (H Q = 72 kms -1 Mpc -1 , SIm = 0.3, 
£Ia — 0.7). This epoch corresponds well to the re-ionization era of the universe. 



Subject headings: galaxies: 
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Introduction 



quasars — galaxies: high-redshift 



laxies: Fell emission — galaxies: 



There is growing evidence that galaxy forma- 
tion and intense star formation are closely re- 
lated to quasar activity, particularly in the early 
universe. For example, the detections of large 
amounts of molecular gas and dust (M mo ; ec . ~ 
1O 1O M , M dust ~ 5 x 1O 8 M ) around quasars 
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at redshifts z > 4 (Isaak et al. 1994; Ohta et 
al. 1996; Andreani et al. 1999; Carilli et al. 2000, 
2002; Omont et al. 1996, 2001) establish a close re- 
lationship between quasars and massive galaxies. 
There is also solid evidence that the co-moving 
space density of quasars and the cosmic star for- 
mation rate both show a strong increase toward 
higher redshifts (Lilly et al. 1996; Connolly et al. 
1997; Tresse & Maddox 1998; Steidel et al. 1999; 
Lanzetta et al. 2002). At low redshifts it has 
been found that characteristic physical properties 
of quasar host galaxies, e.g., bulge mass, luminos- 
ity and stellar velocity dispersion, are well cor- 
related with the mass of the central black hole 
(Kormendy & Richstone 1995; Gebhardt et al. 
2000; Merritt & Ferrarese 2001; Ferrarese & Mer- 
ritt 2001; Tremaine et al. 2002). Moreover, it has 
become clear that all massive galaxies today har- 
bor super massive black holes (SMBHs) in their 
centers, implying that at one time (early in their 
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evolution) these galaxies were hosts to active, i.e. 
quasar or quasar-like, nuclei. 

Since quasars are among the most luminous 
objects in the universe they can be observed at 
practically any distance at which these objects 
occur. This practical aspect together with their 
close relation to galaxy and star formation history 
makes quasars excellent probes to study these pro- 
cesses in the early universe. Assuming that the gas 
of the broad emission line region (BELR) origi- 
nates from the interstellar medium of the quasar 
host galaxies, these lines provide valuable infor- 
mation on the chemical composition and enrich- 
ment history of the quasar host galaxies. Analysis 
of ultraviolet broad emission lines (BELs), as well 
as intrinsic associated absorption lines in quasar 
spectra, indicate that the gas-phase metallicities 
near quasars are at least solar to several times so- 
lar at even the highest measured redshifts (e.g., 
Hamann & Fcrland 1992, 1993, 1999; Ferland et 
al. 1996; Korista et al. 1996; Hamann 1997; Diet- 
rich ct al. 1999; Pettini 1999; Dietrich & Wilhclm- 
Erkens 2000; Warner ct al. 2002; Dietrich ct al. 
2003a, b). These results show that well before the 
epoch corresponding to z ~ 3 to 4 significant star 
formation must have taken place in the galactic 
or proto-galactic cores where these quasars reside 
(e.g., Cen & Ostriker 1999) to enrich the gas phase. 

The epoch of first intense star formation activ- 
ity in quasar host galaxies can also be dated ex- 
ploring the abundance of a-process elements rel- 
ative to iron. According to present chemical en- 
richment scenarios, a-element nuclei are produced 
predominantly in type II SNc with massive pro- 
genitors on time scales of T evo i ~ 2 — 10 Myr. 
On the other hand, the dominant source of iron 
is assumed to be type la SNe at the end point of 
the evolution of intermediate mass stars in binary 
systems, about t pvo i ~ 1 Gyr after the onset of 
the star formation epoch (e.g., Tinsley 1979; Mat- 
teucci & Greggio 1986; Wheeler et al. 1989; Yoshii 
et al. 1996). The amount of iron returned to the 
interstellar medium in SNII ejecta is rather low 
(Yoshii et al. 1996, 1998). It is generally assumed 
that the significantly different time scales of the 
release of a-elements and iron to the interstellar 
medium results in a time delay of the order of 
~ 1 Gyr in the iron enrichment. However, recent 
studies by Matteucci (1994), Friaga & Terlevich 
(1998), and Matteucci & Recchi (2001) indicate 



that this delay may be ~ 0.2 — 0.6 Gyr for ellipti- 
cal galaxies. Detecting Fell emission at high red- 
shift comparable to the relative strength observed 
in quasars at lower redshift indicates that the for- 
mation of the stars which had released the iron 
had occurred ~ 0.3 — 0.8 Gyr earlier in the quasar 
host galaxy. Therefore, the line ratio of a-elcmcnt 
vs. iron emission can be used as a cosmological 
clock (e.g., Hamann & Ferland 1999; Matteucci & 
Recchi 2001). 

A suitable spectral feature for representing the 
presence of a-elements is the prominent Mgn 
resonance doublet at AA2795,2803A (hereafter 
MgnA2798). The generally strong UV Fen multi- 
plets around Mgn in quasar spectra (AA ~ 2000 — 
3000 A; Wills et al. 1980; Grandi 1981; Netzer et 
al. 1985), are probably the most promising for 
obtaining at least some information on the rela- 
tive iron abundance (e.g., Wamplcr & Oke 1967; 
Wills et al. 1985; Hamann & Ferland 1999). Since 
Mg + and Fe + have similar ionization potentials 
(7.6 eV and 7.9 eV, respectively), it is generally as- 
sumed that Mg nA2798 and the Fe n emission both 
originate in the partially ionized zone of the line 
emitting gas. However, it is important to keep 
in mind that the MgnA2798 and Fen emission 
are differently affected by the conditions of the 
line emitting gas. As Wills et al. (1985) pointed 
out, the Fe + zone is much more extended than 
Mg + zone and contains hot regions where magne- 
sium is ionized to Mg ++ , which emits Mgn much 
less efficiently. Furthermore, MgnA2798 and the 
Fe II emission show significantly different radiative 
transfer properties. While MgnA2798 is a pair 
of two resonance lines Fen involves several ten 
thousand lines. Due to these thousands of transi- 
tions, Fe II is subject to processes like line fluores- 
cence, Lya pumping and turbulence (e.g., Verner 
ct al. 1999) which results in a complicated line 
transfer. It is also known that UVFen emission 
suffers destruction by Balmer continuum absorp- 
tion and conversion to optical Fen emission for 
large column densities, Nh ^ 10 24 cm~ 2 , (Wills 
et al. 1985; Joly 1987). Since MgnA2798 forms 
predominantly in regions of higher Balmer con- 
tinuum opacity and these lines are closer to the 
wavelength of the Balmer edge than the strong UV 
Fell emssion, MgnA2798 emission is more reduced 
than UV Fen by Balmer continuum absorption 
(Wills et al. 1985). Another process which affects 
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the Fen and MgnA2798 emission are charge ex- 
change reactions. As pointed out by Joly (1987) 
for T e £j 10 4 K Fe + is ionized to Fe ++ by charge 
exchange reactions while Mg + is nearly unaffected 
by this process. Given the complexity of the 
Fen emission spectrum accurate iron abundances 
are not yet straight forward to deduce and a full 
synthesis of individual quasar spectra is required. 
This seriously complicates such measurements be- 
cause the Fell emission is hard to synthesize (e.g., 
Verner et al. 1999; Collin & Joly 2000; Sigut & 
Pradhan 1998, 2003). 

However, in spite of these uncertainties which 
affect the Fen and the MgnA2798 line emission, 
the Fen/Mgn emission ratio also must reflect the 
Fe/Mg abundance. While, in view of the uncer- 
tainties concerning the detailed physical processes, 
it may not be possible to derive absolute abun- 
dances from these lines, it should, nevertheless, be 
possible to derive at least variations of th Fe/Mg 
ratio with redshift from emission line ratios if ob- 
jects with similar observed physical properties are 
compared. This is feasible owing to the similar 
appearance of their optical-UV spectra across the 
known redshift and luminosity ranges (Dietrich 
et al. 2002b) and the observed similar variability 
properties (Peterson & Wandel 1999, 2000; Onken 
& Peterson 2002), suggesting a generally similar 
physical structure of the AGN emission-line re- 
gions. Assuming that the physical conditions of 
the line emitting gas are comparable for differ- 
ent objects, the effects mentioned above should 
influence the Fen/MgnA2798 line ratio measured 
in individual quasars in a very similar way. To 
minimize the uncertainties that are introduced by, 
for example, different radiative transfer proper- 
ties and different conditions of the line emitting 
gas (T e , NH,n e ), we concentrate (a) on individual 
quasars of similar intrinsic luminosity and (b) on 
composite quasar spectra of comparable luminos- 
ity as the individual high redshift quasars. This 
approach, adopted in this work, reduces the influ- 
ence of a possibly luminosity dependent spectral 
energy distribution (SED) of the incident ioniz- 
ing continuum on the physical conditions of the 
line emitting region. As a consequence, similar 
continua will result in similar conditions of the 
line emitting gas, emphasizing effects of relative 
abundance variations. Assuming that the phys- 
ical parameters affecting Fen/Mgn arc roughly 



the same at all redshifts, then the average of mea- 
sured Fen/Mgn ratios at different redshifts will 
yield useful constraints on the redshift evolution 
of Fe/Mg. 

Observations of the MgnA2798 and of the rest- 
frame UV-multiplets of Fe II in the spectra of sev- 
eral high-redshift quasars have been published by 
Hill et al. (1993), Elston et al. (1994), Kawara et 
al. (1996), Taniguchi et al. (1997), Murayama et 
al. (1998, 1999), Yoshii et al. (1998), Thompson 
et al. (1999), Green et al. (2001), Dietrich et al. 
(2002a), Iwamuro et al. (2002), and Freudling et 
al. (2003). Generally, fairly strong Fen emission 
was found, indicating a very early cosmic epoch of 
the first star formation in the high-redshift quasar 
host galaxies. 

In order to obtain more quantitative informa- 
tion on the beginning of the first star formation 
episode in quasar host galaxies at high redshifts, 
we extend our study of 6 high-redshift quasars 
z ~ 3.4 (Dietrich et al. 2002a) to earlier cos- 
mic epochs. We observed 6 quasars with red- 
shifts 4.4 < z < 5 to extend the redshift range 
and increase the sample size. This high redshift 
range corresponds a time when the universe was 
less than ~ 1.3 Gyr old (H = 72 kms" 1 Mpc -1 , 
n M = 0.3, Q A = 0.7; Carroll et al. 1992; Freedman 
et al. 2001). To put this study on the evolution 
of the Feii/Mgn ratio at high redshifts in a wider 
cosmological context, we analyzed quasar compos- 
ite spectra which we calculated for a wide range in 
luminosity (~ 6 orders of magnitude) and redshift 
(0 i5 z S 5) to study continuum and emission line 
correlations, i.e. the Baldwin-Effect (Dietrich et 
al. 2002b). 

We find a lack of evolution of the measured 
Fell/Mgn emission line ratios in these high- 
redshift quasars compared to those in quasar spec- 
tra at low-redshifts. This indicates that the rel- 
ative iron to magnesium abundance of the gas in 
bright quasars is unchanged since z ~ 5. We con- 
clude that intense star formation activity in the 
host galaxies of the studied high-redshift quasars 
must have started already at an epoch correspond- 
ing to Zf ~ 6 — 9, when the age of the universe 
was ~ 0.5 Gyr. This epoch corresponds well to 
the re-ionization era of the universe. 
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2. Observations and Data Analysis 

We observed 6 quasars with redshifts 4.4 < 
z < 5.0 to measure the Fcn/Mgii ratio in the 
rcstframe ultraviolet (hereafter referred to as the 
z = 4.5 sample). The spectra were recorded using 
infrared spectrographs of the W.M. Keck Observa- 
tory at Hawai'i, the Cerro Tololo Inter- American 
Observatory (CTIO) at Chile, and the ESO VLT 
Paranal Observatory at Chile (Table 1). The 
quasars were selected to be accessible at the time 
of observation and to be bright enough to be ob- 
served within reasonable integration times. Fur- 
thermore, UV emission lines like C IVA1549 are 
uncontaminated by strong and broad absorption 
features. For the selected redshifts the prominent 
Mg nA2798 emission line and most of the strong ul- 
traviolet Fen emission at AA2300 — 2600 A are red- 
shifted into the near infrared H-band. We covered 
practically the whole restframe wavelength range 
AA2200 — 4500A continuously, which made it pos- 
sible to derive reliable continuum fits and, thus, 
at least internally accurate flux values for the ul- 
traviolet Fe II emission. In addition to the z > 4.4 
quasars, Q 2050—358 (z = 3.5) was recorded to en- 
large our sample of quasars at z ~ 3.4 (hereafter 
'the z = 3.4 sample'). The details of the data 
aquisition and processing are described in the fol- 
lowing subsections. 

2.1. Keck II NIRSPEC Observations 

The high redshift quasar PC 1158+4635 (z=4.733; 
Schneider et al. 1989) was observed in service 
mode with NIRSPEC at KECK II in the near 
infrared domain on May 19/20, 2000. The see- 
ing was ~ 0".5, estimated from the spatial pro- 
file of the standard star spectra. The 1024 x 
1024 InSb detector array was used in the longslit 
mode (0".76x 42", 0".19/pixcl). To cover the 
near infrared spectral range three settings of the 
grating were used which yield the wavelength 
ranges (1.26 - 1.54)^m, (1.52 - 1.81) urn, and 
(2.05 — 2.48)/um which cover the traditional J, 
H, and K bands, respectively. Making use of the 
slit length the location of the object along the slit 
was alternated by ~ 20" for subsequent exposures 
to assure a proper sky-correction. Due to the high 
brightness of the night sky in the near infrared, 
the integration times of the individual exposures 
of the quasar were limited to 120 sec and 400 sec, 



respectively. The observation log and the total in- 
tegration times for the individual bands are given 
in Table 1. 

A mean dark frame was calculated using 5 in- 
dividual dark frames. Flatfield frames were taken 
in the lamp-on and lamp-off mode for each wave- 
length range. The frames corrected for ther- 
mal background were used to compute normalized 
flatfield frames for the flatfield correction (Home 
1986). Argon comparison spectra were recorded 
for wavelength calibration for each wavelength 
range. In the J to H infrared wavelength range the 
spectral resolution, measured using the FWHM 
of strong night sky emission lines, amounts to 
R = 570 corresponding to Av ~ 520 kms -1 and 
in the K wavelength range it amounts to R = 650 
(Av ~ 460 kms -1 ). Because the sky brightness 
can vary significantly on time scales of minutes 
(comparable to our exposure times) we used the 
sky spectrum measured along the slit in each frame 
to subtract the sky emission. A 3 rd to 6 th order 
polynomial fit was calculated for each wavelength 
element to describe the spatial intensity distribu- 
tion of the night sky emission. The fit based on 
~ 18" wide regions which were ~ 5" separated 
from quasar or standard star. 

The standard star HD 105601 (m K = 6.685, 
A2V; Elias ct al. 1982; Arnaud et al. 1989) was 
observed for each setting before or after the obser- 
vations of the quasar spectrum for flux calibration 
and correction of the strong atmospheric absorp- 
tion features. To obtain a sensitivity function we 
assumed that the spectral energy distribution of 
HD 105601 can be described in the near infrared 
wavelength range with a black body spectrum. We 
applied T e ff = 8995 K to computed a black body 
spectrum for the observed wavelength range. This 
black body spectrum was scaled to match the ap- 
parent magnitudes of HD 105601 in the J-, H-, and 
K-band (Elias et al. 1982). 

We used the spatial curvature of the standard 
star spectrum to trace the location of the slightly 
curved quasar spectrum. The spectra of the indi- 
vidual exposures were extracted using an optimal 
extraction routine (Home 1986). 

The 1 D quasar spectra were corrected for broad 
atmospheric absorption bands using an appropri- 
ately scaled transmission function derived from 
HD 105601. Cosmic-ray events in the individ- 
ual 1 D spectra were corrected by comparing the 
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spectra with each other. Finally, we computed 
a weighted mean of the quasar spectra for each 
setting, with the weight given by the signal-to- 
noise ratio of the continuum. These mean quasar 
spectra were flux calibrated using the correspond- 
ing sensitivity function for the wavelength range. 
The mean spectra of each setting were rebinned 
to a uniform stcpsize in wavelength and merged. 
The flux calibrated spectrum of PC 1158+4635 is 
shown in Figure 1. 

2.2. CTIO Observations 

On Sept. 14 to 16, 2000 we used the near in- 
frared spectrometer and camera OSIRIS attached 
to the 4 m Blanco telescope at Cerro Tololo Inter- 
American Observatory (CTIO), Chile, to observed 
BRI 0019-1522, BR0103+0032, PSS J0248+1802, 
Q 2050-359, and BRI 2237-0607. The seeing was 
typically ~ 0".8 during these three nights. The 
detector was a 1024 x 1024 Hawai'i HgCdTe ar- 
ray (18.5^m/pixel). We used OSIRIS in the low- 
resolution cross-dispersed mode with a longslit of 
1".2 x 30". The pixel scale was 0".403/pixel. In 
the cross-dispersed mode the J, H, K bands are 
recorded in adjacent orders in a single exposure 
yielding the entire near infrared wavelength range, 
i.e., (~ 0.9 - 2.5)/im. Making use of the long-slit 
mode we placed the object at seven different loca- 
tions along the slit, separated by 3" each for sub- 
sequent exposures to optimize the sky-correction. 
For flux calibration and correction of the strong 
atmospheric absorption bands, separating the J- 
and H-band and the H- and K-band, respectively, 
and additional weaker atmospheric (absorption) 
features we observed standard stars several times 
during each night (Elias ct al. 1982; Perryman et 
al. 1998). In Table 1 the exposure times for the 
quasars are given. 

For each of the different individual exposure 
times of the science frames, i.e., quasars (200 sec, 
300 sec) and standard stars (3 sec to 20 sec), we 
recorded 20 to 25 dark frames. They were used to 
calculate average dark frames for each specific ex- 
posure time. The quasar and standard star spec- 
tra were dark corrected, subtracting the average 
dark frame of the corresponding exposure time. 

Flatficld frames were taken in the lamp-on and 
lamp-off mode, both with 200 and 300 seconds in- 
tegration time, respectively. The average lamp- 
off mode frames were subtracted from the aver- 



age lamp-on mode frames to correct for the ther- 
mal background. These thermal background cor- 
rected mean flatfield frames were split up into sep- 
arate 2D spectra for the J, H, and K band spec- 
trum. Those were used to compute normalized 
flatfield frames (Home 1986) which we employed 
for 2 D flatficld correction for each individual spec- 
tral band. 

For wavelength calibration we took helium- 
argon-ncon comparison spectra. The location of 
the emission lines in the individual calibration 
frames differ by less than £ 0.1 pixel compared 
to each other. Therefore, an average calibration 
frame was computed and was used for the wave- 
length calibration. The uncertainty of the wave- 
length calibration is AA ~ 0.3 A. The spectral 
overlap of the individual spectral bands is of the 
order of ~ 350 A in the quasars restframe which 
allow a solid intercalibration to obtain a contin- 
uous spectrum of the near infrared wavelength 
range. The achieved spectral resolution amounts 
to R ~ 1300 in J, H, and K-band, respectively, 
based on the measured FWHM of strong night sky 
emission lines and isolated lines of the wavelength 
calibration frames. 

We corrected each frame individually for the 
night sky intensity. For this purpose a 3 rd order 
polynomial fit was calculated for each wavelength 
element to fit the spatial intensity distribution of 
the night sky emission. The fit was based on two 
regions, 2" to 20" wide, depending on the location 
of the object along the slit. These regions were 
separated by <~ 10". 

To derive a sensitivity function we assumed that 
the IR spectral energy distribution of a A-type 
and G-type star can be described with a black 
body spectrum of a specific T e //- We used the 
black body temperatures given by Kurucz (1992) 
for the different spectral types of the stars we 
observed. We applied T eff = 5850 K (HD 1205, 
HD 25402, HD 198678), T eff = 8250 K (HD 19904, 
HD 205772), and T eff = 8790K (HD2811) to 
compute a black body spectrum for the observed 
wavelength range. The black body spectra were 
scaled to match the apparent magnitudes of the 
observed standard stars in the J-,H-, and K- 
band. For HD2811, HD 19904, and HD 205772 
we used the J, H, and K magnitudes given by 
Elias et al. (1982), while for HD 1205, HD 25402, 
and HD 198678 we had to calculate the apparent 
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near infrared magnitudes based on the apparent 
V-band magnitude. A sensitivity curve was calcu- 
lated for each standard star. In general, the sen- 
sitivity functions given by the individual stars are 
nearly identical in shape and strength, within less 
than 2% . However, for a ~ 200 A wide region at 
the beginning and the end of each spectral wave- 
length range they differ less than ~ 20% . But 
these regions fall into the strong absorption bands 
which separate the J, H, and K-band. We used the 
individual sensitivity functions obtained from the 
standard stars, we observed, to compute a mean 
sensitivity function for each wavelength range. 

The quasar spectra were extracted using an op- 
timal extraction algorithm (Home 1986). The 
width of the spatial profile for the quasar spec- 
tra was the same as measured for the stars be- 
cause high-redshift quasars can be treated as point 
sources. The individual spectra of the quasars 
were flux calibrated employing the mean sensitiv- 
ity function of the corresponding spectral band. 
To correct for cosmic-ray events the individual 
1 D spectra were compared with each other. A 
weighted mean spectrum was calculated for each 
quasar. The weight was given by the mean signal- 
to-noise ratio in the continuum across the spec- 
trum. The mean quasar spectra were corrected for 
atmospheric absorption using appropriately scaled 
transmission functions provided by observed spec- 
tra of the standard stars. The flux calibrated 
quasar spectra observed with OSIRIS at the 4 m 
Blanco telescope are shown in Figure 1. 

2.3. Paranal VLT/ISAAC Observation 

The high-redshift quasar SDSS 1204-0021 (z = 
5.03; Fan et a. 2000) was observed in service mode 
with ISAAC attached to the VLT UT1 Antu in the 
H-band spectral range. The seeing was estimated 
from the spatial profile of standard star spectra 
to - l".l (April 11, 2001) and - 1".2 (June 14, 
2001). A Hawai'i 1024 x 1024 array from Rock- 
well was used in the longslit mode (1" x 120", 
0".147/pixel). We varied the location of the ob- 
ject along the slit by 60" (with an additional ran- 
dom offset within 5" for subsequent exposures to 
optimize the sky correction). For flux calibration 
and correction of the strong atmospheric absorp- 
tion bands at the beginning and end of the H-band 
spectrum a standard star was observed each night. 
On April 10 and June 14, 2001, respectively, ten 



quasar spectra with 180 sec exposure time were 
recorded which results in a total exposure time of 
1 hr (Table 1). 

For each night three dark frames with 5 sec 
(mean of 10 individual exposures) and 180 sec ex- 
posure time (mean of 3 individual exposures) were 
recorded. These frames were averaged for the cor- 
responding exposure time and used for dark cor- 
rection. Flatfield frames were taken in the lamp- 
on and lamp-off mode for each night. The lamp- 
off mode frames were subtracted from the adja- 
cent lamp-on mode frame. The resulting flatfield 
frames were averaged and were normalized to ap- 
ply a 2 D flatfield correction. 

Argon-xenon comparison spectra were recorded 
for wavelength calibration. The wavelength range 
amounts to ~ 14175 to 19040 A (April 10) and 
~ 14195 to 19065 A (June 14) with a stepsize of 
<~ 4.8 ± 0.5 A/pixel. The wavelength calibration is 
based on 20 individual lines. The spectral resolu- 
tion, measured using the FWHM of strong night 
sky emission lines and isolated lines in the argon- 
xenon spectra, amounts to R ~ 570. 

We subtracted the night sky intensity for each 
spectrum individually. A 3 rd order polynomial 
fit was calculated for each wavelength element to 
describe the spatial intensity distribution of the 
night sky emission. The fit based on two regions, 
~ 25" and ~ 80" wide, which were separated by 
~8". 

To derive a sensitivity function the standard 
stars HD 87097 (April 10) and HD 143567 (June 
14) were observed. We used the black body tem- 
peratures given by Kurucz (1992) with T e ff = 
5850 K (HD 87097, G3V) and T eff = 10700 K 
(HD 143567, B9V) to compute a black body spec- 
trum for the H-band range. These black body 
spectra were scaled to match the apparent magni- 
tudes of the observed standard stars. A sensitivity 
curve was calculated for each standard star based 
on the V-band magnitude. 

The spectra of the quasars which can be re- 
garded as point sources, were extracted using an 
optimal extraction algorithm. The width of the 
spatial profile for the quasar spectra was the same 
as measured for the standard stars. 

The individual spectra of SDSS 1204-0021 
were flux calibrated employing the correspond- 
ing sensitivity function. To correct for cosmic-ray 
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events the individual 1 D spectra were compared 
with each other. For each quasar a weighted mean 
spectrum was calculated. The weight was given 
by the mean signal-to-noise ratio in the continuum 
across the spectrum. The mean quasar spectra 
were corrected for atmospheric absorption using 
appropriately scaled transmission functions pro- 
vided by observed spectra of the standard stars. 
Finally, the spectra obtained in April and June 
2001 were averaged according to the signal-to- 
noisc ratio in the continuum. In Figure 2 the flux 
calibrated H-band spectrum of SDSS 1204-0021 
is shown. 

2.4. Composite Quasar Spectra 

The analysis of composite spectra has the ad- 
vantage of more clearly representing the average 
properties of the sample used to calculate each 
composite. We recently compiled a large sample of 
rest-frame ultraviolet and optical spectra for more 
than 800 type 1 AGNs (Dietrich et al. 2002b). 
About ~ 300 of these spectra cover the ultravi- 
olet Fen - Mgn range (AA2000 - 3000 A). A large 
fraction of the spectra were obtained by several 
groups for different studies over the last 20 years 
using ground-based instruments as well as Inter- 
national Ultraviolet Explorer (IUE) and Hubble 
Space Telescope (HST). In addition, we retrieved 
all publicly available AGN spectra from the HST 
MAST archive. The type 1 AGNs (i.e., Seyfert 1 
and quasars) of the sample cover a redshift range 
of < z < 5 and nearly 6 orders of magnitude in 
luminosity. We excluded Broad- Absorption Line 
quasars (BALQSOs) when calculating the com- 
posite spectra, although there are indications that 
BAL QSO emission line properties do not differ 
from non-BAL quasars (Weymann et al. 1991). 

We analyzed these composite spectra which we 
calculated to investigate continuum vs. emis- 
sion line relations (Dietrich et al. 2002b) for the 
Feil/Mgil ratio. We selected a narrow luminosity 
range to calculate composite spectra for different 
redshift bins with nearly constant luminosity (the 
redshift sample). The location of this luminosity 
range was chosen to achieve a comparable number 
of individual quasars contributing to each redshift 
bin in the z - AL A (1450A) plane. We used a lumi- 
nosity range of 46.16 < log AL A (1450A) < 47.16 
ergs" 1 (i.e., 43.0 < log L A (1450A) < 44.0). This 
range contains 323 quasars with z > 0.5 and an av- 



erage luminosity of log AL A (1450A) = 46.67±0.12 
ergs -1 . We computed composite spectra for red- 
shift bins with Az = 0.5 centered on redshift 
z = 0.75, 1.25, ...,4.25, and4. 75, respectively (Di- 
etrich et al. 2002b). The composite spectra char- 
acterized by nearly constant luminosity provide 
valuable Fen/Mgn measurements as a function 
of redshift. The quasar sample was also divided 
into 11 subsets with respect to the luminosity 
-L A (1450A) (the luminosity sample, see Table 2). 
Except for the lowest luminosity bin, each com- 
posite spectrum covers a range of A log XL\ = 0.5 
in luminosity starting at log AL A (1450A) = 43.16 
ergs -1 . The composite spectra with increasing lu- 
minosity yield information on the Feil/Mgil ratio 
as a function of luminosity. 

3. Fitting of the Quasar Spectra 

To measure the Fe n and Mg n emission strengths, 
we first transformed the individual high-rcdshift 
quasar spectra into the restframe. Although the 
MgnA2798 emission line doublet sometimes ap- 
pears to be relatively isolated in quasar spectra, 
the superposition of many thousands of discrete 
ultraviolet Fell emission lines forms broad emis- 
sion blends (e.g., Wills et al. 1985, hereafter 
WNW85; Verner et al. 1999; Sigut & Pradhan 
1998, 2003), which tend to blend with the broad 
base of Mgll complicating the measurement of 
Mgn and Fen. Because of the line blending and 
the related difficulty defining the continuum level, 
these lines cannot be measured individually. How- 
ever, as demonstrated by WNW85, it is possible 
to obtain the Fell emission strength in a multi- 
component fit approach. Therefore, we assumed 
our observed quasar spectra (z=3.4 and z=4.5 
sample), as well as the composite quasar spec- 
tra to consist of a superposition of the follow- 
ing four components, (i) a power law continuum 
(F„ ~ v a ), (ii) Balmcr continuum emission, (iii) 
a pseudo-continuum due to merging Fe II emission 
blends, and (iv) an emission spectrum of other in- 
dividual broad emission lines, such as Mgll. The 
power law continuum fit to the individual and 
composite quasar spectra was iteratively deter- 
mined taking into account the measured strength 
of the Fe II emission and Balmcr continuum emis- 
sion derived by the fitting process. We already 
applied this multi-component fit approach suc- 
cessfully to a small sample of high redshift quasars 



7 



at z ~ 3.4 (Dietrich et al. 2002a). 

3.1. Non-stellar Continuum 

We first determined the underlying non-stellar 
power law continuum, component (i), from spec- 
tral windows which arc free (or almost free) of con- 
tributions by the components (ii) to (iv) . For this 
purpose it is very important that the spectrum 
covers a wide wavelength range. The high-redshift 
quasars we observed allow us to use a 20 A wide 
continuum window, centered at A c ~ 4000 ± 50 A 
which is nearly free of emission line contamination 
and hydrogen Balmer continua. Only minor Fe II- 
cmission can be expected for this wavelength re- 
gion (WNW85; Verner et al. 1999; Corbin & Boro- 
son 1996) and hence will not significantly affect the 
continuum setting. For 4 of the individual high- 
redshift quasars (BRI 0019-1522, BR 0103+0032, 
PC 1158+4635, BRI 2237-0607) we have also rest- 
frame UV spectra. These spectra provide con- 
tinuum windows at AA ~ 1330 — 1380 A, and 
AA ~ 1440-1470 A, and AA ~ 1685-1695 A which 
are nearly uncontaminatcd by line emission as well 
(Francis et al. 1991). For those quasars for which 
we had no access to the short wavelength UV part 
of the spectrum (PSS J0248+1802, Q 2050-358) 
we used continuum windows at ~ 2230 A to esti- 
mate the continuum strength. The UV data al- 
low us even better constraints on the non-stellar 
continuum strength for these objects. The un- 
certainty introduced by estimating the continuum 
level, based on at least two of the spectral ranges 
described above, is estimated to be of the order of 
less than - 10% . 

3.2. Balmer Continuum 

To estimate the Balmer continuum emission 
spectrum for our fits we assumed gas clouds of uni- 
form temperature (T e = 15000 K) which are par- 
tially optically thick. In this case the Balmer con- 
tinuum spectrum (A < 3646 A) can be described 

by 

F BaC = F BE Bx{Te){l _ e - TA); A < 

with B\(T e ) as the Planck function at the elec- 
tron temperature T e , t\ as optical depth at A, 
and F BE as a normalized estimate for the Balmer 
continuum flux density at the Balmer edge at 



A = 3646 A (Grandi 1982). The strength of the 
Balmer continuum emission can be estimated from 
the flux density at A ~ 3675 A, after subtraction 
of the power-law continuum component, since at 
this wavelength there is no significant contami- 
nation by Fen emission (WNW85; Verner et al. 

1999) . The A3675A restframe flux density level 
was therefore used to normalize the Balmer con- 
tinuum spectrum. At wavelengths A > 3646 A 
higher order Balmer lines arc merging to a pseudo- 
continuum, yielding a smooth rise to the Balmer 
edge (WNW85) . We used the results of the model 
calculations provided by Storey & Hummer (1995) 
(case B, T e = 15000K, n e = 10 8 - 10 10 c m - 3 ) to 
estimate the strength of high order Balmer emis- 
sion lines, assuming a Gaussian profile for each 
high order Balmer line (FWHM = 3000 kms" 1 ). 
We calculated several Balmer continuum spectra 
for T e = 15000 K and 0.1 < r A < 2 to ob- 
tain Balmer continuum emission templates. These 
Balmer continuum templates were supplemented 
for A > 3646 A with high order Balmer emission 
lines with 10 < n < 50, i.e., H$ and higher. 

3.3. Fen Emission 

Calculating the Fe n emission spectrum is much 
more difficult and the influence of unknown pa- 
rameters such as internal turbulence velocities, 
emission line transport, pumping by the incident 
continua, line fluorescence, and metallicity, which 
affect the emergent Fen spectrum, are still not 
well understood (e.g., WNW85; Netzer et al. 1985; 
Joly 1987; Bautista & Pradhan 1998; Sigut & 
Pradhan 1998; Verner et al. 1999; Collin & Joly 

2000) . In spite of these uncertainties it has been 
shown that the Fe n emission strength of Seyfert 1 
galaxies and quasars can be measured using em- 
pirically derived Fen emission template spectra 
(WNW85; Laor et al. 1997; Mcintosh et al. 1999; 
Dietrich et al. 2002a; Iwamuro et al. 2002). We 
fitted the Fe n emission in our quasar spectra us- 
ing scaled and broadened empirical Fe n emission 
template spectra to derive relative Fell emission 
strength values. For the ultraviolet wavelength 
range these templates had been carefully extracted 
from HST observations of I Zwl by Vestergaard & 
Wilkes (2001). These emission templates include 
contributions of Fe n and Fe III. In the fitting pro- 
cess the strength and the broadening of the Fe 
emission template are free parameters. 



8 



The use of different empirical Fe n emission 
templates, in particular in the range of the broad 
MgnA2798 emission line, introduces additional 
uncertainty to measured Fe n/Mg n ratios. For ex- 
ample, the empirical Fe emission template we use 
in our analysis underestimates the Fe flux in the 
region. Vcstcrgaard & Wilkes (2001) set the flux 
level beneath the broad MgnA2798 emission line 
to zero in the Fe template extracted from IZwl. 
This approach is reasonable since detailed model 
calculations (Verner et al. 1999; Sigut & Pradhan 
2003) indicate that the Fe flux level in this wave- 
length range is of the order of only 15 to 20 % 
compared to neighboring Fe emission. The treat- 
ment of the Fe emission beneath Mgn line has a 
major impact on the Mgn flux. Since the broad 
component can contain <~ 50—80 % of the total line 
flux, by neglecting the presence of the broad pro- 
file component can easily cause an underestimate 
of the Mg II flux by a factor of ~ 2 (Vestergaard 
& Wilkes 2001; Dietrich et al. 2002b). 

We explored the effects of non-zero Fe n fluxes 
beneath the Mgn emission line. The gap in the 
Fe emission template, we used for the analysis 
(Vestergaard & Wilkes 2001), was replaced by 
a low order polynomial (2 nd order) which mini- 
mal flux level accounts to ~ 15 % of the imme- 
diately close by located Fe emission peaks. This 
choice is motivated by the photoionization calcula- 
tions mentioned above (Verner et al. 1999; Sigut & 
Pradhan 2003). This additional flux increases the 
integrated Fen flux in the range AA2200 - 3090 A 
by less than 2%. Hence, the integrated Fen flux 
is nearly unchanged by filling the gap beneath the 
Mgn line with a more realistic estimate of the 
complex Fe II emission. 

To estimate the implications of the modified 
Fe emission template on the measured Mg nA2798 
emission line strength we re-analyzed the Mg n line 
profile in each individual high-redshift and com- 
posite quasar spectra. The line profile was fit with 
a narrow and a broad Gaussian component. The 
narrow component was nearly unchanged. The 
strength of the broad component was reduced by 
~ 6 % on average. The cumulative effect of the 
slightly increased Fen emission and the reduced 
strength of the Mgn emission flux amounts to 
~ 10% increase for the measured Feil/Mgn ra- 
tios, nearly independent of the Mgn profile width. 



3.4. Strong Broad Emission Lines 

The broad emission line profile of Mg nA2798 
was fit after subtracting the power law continuum, 
the Fell emission, and the Balmer continuum 
emission. Generally, the MgnA2798 emission line 
profile could be reconstructed with two Gaussian 
components, a narrow component (FWHM(Mgn) 
= 2270±570 kms -1 ) and a blueshifted broad com- 
ponent (Av ~ -1000 kins" 1 , FWHM(Mgn) = 
6300± 1300 kms -1 ). The approach to reconstruct 
the Mgn emission line profile with two Gaussian 
components was only chosen to measure the line 
flux; each individual component has no physical 
meaning by itself. 

3.5. Internal Reddening 

In recent years growing evidence has ap- 
peared for the presence of large amounts of dust 
(Mdust ^ 1O 8 M ) in the host galaxies of high- 
redshift quasars (e.g., Guilloteau et al. 1997, 1999; 
Carilli et al. 2000; Omont et al. 2001). It is as- 
sumed that the dust is distributed in a kiloparsec- 
scale warped disk (Sanders et al. 1989) which is 
illuminated by the central AGN. The observed 
dust emission spectra from 3 to 30 /Ltm can be ex- 
plained by such a model as shown by Andreani et 
al. (1999) and Willott et al. (2000). However, since 
the presented spectra are typical quasar spectra 
with prominent broad emission lines, it is unlikely 
that the broad-line region (BLR) is significantly 
blocked by dust (e.g., Netzer & Laor 1993). Even 
if the radiation has to pass through an (external) 
dust screen the extinction of the Mgn and the 
UV Fe n-emission (having about the same mean 
wavelength) will be comparable and the line ratio 
is not expected to be significantly modified. 

3.6. Multi-Component Fit of the Quasar 
Spectra 

We have analyzed the individual high-redshift 
quasar spectra, as well as the composite quasar 
spectra in a multi-component fit approach us- 
ing emission templates as described above. The 
spectra are transformed to their restframe using 
the corresponding redshift given in Table 2. The 
quasar spectra were iteratively fitted first with a 
power law continuum, then a Balmer continuum 
emission template and an Fe emission template 
were employed, and finally a two Gaussian pro- 
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file fit was used to measure the broad Mg nA2798 
emission line. The strength of the Fe emission 
templates was varied, as well as the width of the 
individual iron features. The broadening of the 
empirical Fe emission template (Vestergaard & 
Wilkes 2001), was achieved by convolution with 
a Gaussian profile. The width of the iron emission 
features of the best fit is consistent with the pro- 
file width of the CivA1549 and MgnA2798 emis- 
sion lines to within ~ 1000 kms -1 . We also var- 
ied the strength and the optical depth t\ of the 
Balmer continuum emission (0.1 < t\ < 1.0). 
Minimal x 2 fits were calculated to determine the 
best fit. In order to estimate the accuracy of our 
fits, we also calculated fits for settings of differ- 
ent parameters around the best fit values. From 
the widths of the resulting x 2 distributions we es- 
timated mean errors for these parameters which 
were employed to estimate the uncertainties of the 
flux measurements. In spite of a signal-to-noise 
ratio of S/N ~ 10 in the continuum of the indi- 
vidual high-redshift quasar spectra, the UV Fell 
and Mgn emission strength can be measured to 
within - 12 - 20% (Fen) and - 10% (Mgn), re- 
spectively. The errors of the Fen/Mgn ratio de- 
termined for the composite spectra also take into 
account the uncertainties introduced by the indi- 
vidual spectra that contribute to the correspond- 
ing composite spectrum. We estimated this uncer- 
tainty from the rms-spectrum which is computed 
for each composite spectrum. 

4. Results 

In Figure 1 we present the quasar spectra of the 
z = 4.5 sample together with the power law con- 
tinuum fit, the appropriately scaled Fe and Balmer 
continuum emission templates, and the fit of the 
MgllA2798 emission line. The achieved multi- 
component fit results benefit significantly from the 
large covered wavelength range to define a reli- 
able continuum fit, in contrast to prior studies. 
The residuum, subtracting the best fit from the 
quasar spectrum, is displayed in the bottom panel 
for each quasar. The remaining quite strong fea- 
ture in the residuum spectra at A ~ 3200 A is due 
to prominent iron emission blends, multiplcts M6 
and M7, that is not covered by the Fe template 
which we are using. The ultraviolet Fe n emission 
line flux was measured from the scaled Fe tem- 
plates for a wavelength range of AA2200-3090A. 



The integrated flux of MgnA2798 was determined 
from the two component Gaussian fit. The mea- 
sured flux for both emission features are given in 
Table 2. 

For SDSS 1204-0021 only the H-band spec- 
tral range was observed. Instead of applying an 
uncertain multi-component fit - the continuum 
strength and the contribution of the Balmer con- 
tinuum emission can not be reliably determined - 
we compared the strength of the MgnA2798 and 
the broad Fen emission feature at A ~ 2500 A 
with the quasars at z ~ 4.5. We calculated a 
mean z ~ 4.5 quasar spectrum based on the other 
quasars of the z = 4.5 sample. In Figure 2 the 
spectrum of SDSS 1204-0021 transformed to the 
rcstframe is shown together with the mean z ~ 4.5 
quasar spectrum. This mean high-redshift quasar 
spectrum was scaled by multiplying with a con- 
stant to match the spectrum of SDSS 1204-0021. 
The difference spectrum is displayed in the bot- 
tom panel of Figure 2. The overall spectral slopes 
of both these spectra are very similar. Since the 
spectra were scaled to match the overall spec- 
tral shape a small residuum is left of Mg n emis- 
sion. The Mg II emission line profile is significantly 
broader in SDSS 1204-0021 (FWHM ~ 7400 
kms -1 ) than in the mean z ~ 4.5 quasar spectrum 
(FWHM ~ 3000 kins" 1 ) and it does not show 
a pronounced narrow component. Although the 
profiles of Mg n look quite different this is predom- 
inantly due to the missing narrow component in 
SDSS 1204—0021, as can be seen in the difference 
spectrum (Fig. 2, bottom panel). The residual flux 
in the difference spectrum amounts to ~ 5 % of 
the total Mgn flux in both quasar spectra. Al- 
though the Feil/Mgn ratio can not be measured 
in a multi-component fit, the comparison with the 
mean z ~ 4.5 quasar spectrum indicates that the 
ratio should be very similar to the mean value of 
the individual quasars at z > 4. 

4.1. Fen/Mgn Ratio as a Function of Red- 
shift 

First, we compare the Fen/Mgn ratios for the 
z = 3.4 and z — 4.5 quasar sample we observed 
for individual quasars. The measured Feli/Mgn 
ratios which are given in Table 2, are displayed 
as function of redshift in Figure 3. The results 
for the quasars of the z — 3.4 sample (Dietrich 
ct al. 2002a) are supplemented with Q 2050— 359 
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(z = 3.51) which we observed at CTIO. There 
is obviously no strong indication for an evolution 
of the Fen/Mgll ratio from z ~ 3.2 up to z ~ 
4.8. Including the result of the comparison of the 
SDSS 1204-0021 spectrum with the mean z = 4.5 
quasar spectrum, the Feil/Mgll ratio shows no 
variation up to z ~ 5. Although our observation 
of PC 1158+4635 indicates a lower Feil/Mgll ra- 
tio, SDSS 1204-0021 and recent NICMOS/HST 
observations of three z ~ 6 quasars show no de- 
clining Feil/Mgll ratio (Freudling et al. 2003). 

To investigate the evolution of the Feil/Mgll 
ratio over a wide range in cosmic time, particu- 
larly from the local universe to z ~ 5, we measured 
this emission line ratio in the quasar spectra of the 
redshift sample. The criterion of nearly constant 
and also comparable luminosity of the composite 
spectra of this sample (Table 2) compared to the 
z = 3.4 and z — 4.5 quasar samples minimizes lu- 
minosity effects on the emission line strength like 
the Baldwin effect (e.g., Osmcr & Shields 1999; 
Croom et al. 2002; Dietrich et al. 2002b; Shang et 
al. 2003). The MgnA2798 and UV Fen line fluxes 
we measure for the quasars of the redshift sample 
are given in Table 2. In Figure 3 these Fen/Mgn 
ratios are plotted together with the results for 
the individual high- redshift quasars. The analy- 
sis of the redshift quasar spectra provides further 
evidence that the ratio of iron to a-elements ex- 
pressed by Feil/Mgll lacks an evolution for lumi- 
nous quasars from the local universe (z ~ 0) up to 
z ~ 5, i.e., an cosmic age of ~ 1 Gyr (Fig. 3). 

4.2. Feii/Mgn Ratio as a Function of Lu- 
minosity 

To investigate the dependence of the Fe n/Mg n 
ratio on luminosity we analyzed the composite 
quasar spectra of the luminosity sample which 
cover nearly <~ 6 orders of magnitude in lumi- 
nosity (Dietrich et al. 2002b). In Figure 4 the 
Feil/Mgll ratios measured for these composite 
spectra (luminosity sample) are displayed as a 
function of continuum luminosity ALa(1450). In 
addition, the Feil/Mgll ratios of the individual 
high-redshift quasars of the z = 3.4 and z = 4.5 
samples are shown. Within the errors there is a 
weak trend for a higher Feil/Mgll ratio towards 
higher luminosities (Fig. 4) . It appears that the 
ratio is nearly constant for logX La(1450A) < 45 
ergs -1 with Fell/Mgll= 2.1 ±0.3 and amounts to 



Fell/Mgll= 4.0 ± 0.8 at higher luminosities. The 
composite quasar spectrum at log\L\ = 43.51 
ergs -1 consists of only 3 individual spectra from 
which two cover the UV Fen - Mgn wavelength 
range. One of these two quasars is WPVS007, a 
well-known narrow-line Seyfcrt 1 (NLS1) galaxy. 
The HST spectrum that we retrieved from the 
HST data archive, shows exceptionally strong Fell 
emission (Grupe et al. 1995) which is a common 
characteristic of NLS1 galaxies (e.g., Mathur 2000; 
Pogge 2000). Thus WPVS007 dominates the com- 
posite spectrum at log\L\ = 43.51 ergs -1 and 
results to the high value of Fen/Mgn = 5.6 ± 1.2. 

4.3. Comparison with other Studies 

During the last few years there have been a 
few similar studies of the UV Feil/Mgll ratio 
in quasars at high rcdshifts (e.g., Kawara et al. 
1996; Taniguchi et al. 1997; Murayama et al. 
1999; Thompson et al. 1999; Iwamuro et al. 2002; 
Freudling et al. 2003), as well as for lower redshift 
samples (Wills et al. 1980; Grandi 1981; Wills et 
al. 1985). However, to compare those results with 
the Feil/Mgll ratios we measure here (z = 3.4, 
z = 4.5, and redshift and luminosity quasar sam- 
ples), it is necessary to adjust the published ra- 
tios to reflect the same spectral range we used 
to measure the Fen emission (see below). The 
comparison is further complicated by the differ- 
ent Fen emission templates which were used to 
measure the Fen flux, as well as taking into ac- 
count the contribution of the Balmcr continuum. 
Wills et al. (1985) applied an Fen emission tem- 
plate based on their photoionization model calcu- 
lations (Netzer & Wills 1983), while Thompson 
et al. (1999) and Iwamuro et al. (2002) used the 
LBQS mean quasar spectrum (Francis et al. 1991) 
to extract an empirical Fe n emission template. 
Another more serious uncertainty is introduced 
by the generally restricted wavelength coverage of 
prior studies that made a proper definition of the 
underlying continuum strength very difficult, in 
particular, in the restframe UV wavelength range 
of high-redshift quasars. The quasar spectra Wills 
et al. (1985) analyzed covered at least the spec- 
tral range from Cm]A1909 up to H/3A4861. The 
high redshift quasars (z = 3.4 and 4.5 samples) in 
the Thompson et al. (1999) study cover ~ 600 A 
around the MgnA2798 emission line and, hence, 
does not cover the entire Fen bump. The in- 
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dividual spectra of the Thompson et al. sample 
have been re-analyzed by Iwamuro et al. (2002), 
who also observed the Mg n - UV Fe n wavelength 
range for 13 quasars with z > 4.4, across a wave- 
length range of <~ 1200 A. Iwamuro et al. (2002) 
also presented measurements of Feil/Mgn for the 
early data release of the SDSS QSO survey, as well 
as for quasar spectra retrieved from the HST and 
IUE data archives for a local quasar population. 
Recently, Freudling et al. (2003) presented results 
of measured Fen/Mgil ratios for three quasars at 
z ~ 6. Although the uncertainties are quite large 
the Fe ii/Mgn ratios are comparable to the results 
of the z > 4 quasars both here and by Iwamuro et 
al. (2002). 

To reduce the most obvious source of differ- 
ences for the individual published Feil/Mgn ra- 
tios we adjusted the Fen/Mgn emission line ra- 
tios of prior studies to the integration window 
that we used to measure the Fell emission line 
flux. Therefore, we compared the Fe emission flux 
we measured for the integration range we used 
(AA2200 - 3090 A) with the those resulting by ap- 
plying different integration windows. (WNW85; 
Thompson et al. 1999; Iwamuro et al. 2002). 

In Figure 5 we present our measurements of the 
Fell/Mgn ratio for the z = 3.4, z = 4.5, and red- 
shift samples in comparison with the rescaled re- 
sults presented by Wills et al. (1985), Thompson 
et al. (1999), Iwamuro et al. (2002), and Freudling 
et al. (2003). At higher redshifts (z > 3) the in- 
tercalibrated Fe ii/Mgn ratios are quite consistent 
although the scatter of the individual measure- 
ments given by Iwamuro et al. (2002) is large. 
In particular, the Fen/Mgn ratios for the high- 
rcdshift (z > 4) composite spectra of the SDSS 
sample which are also those with the highest lumi- 
nosity of this sample (Iwamuro et al. 2002) agree 
quite well with the results of our composite spectra 
(Fell/Mgn~ 4.4 to 5.7 compared to 5.4±1.3 here, 
Table 2). The Feil/Mgn measurements for three 
quasars at the highest redshift indicate that com- 
parable high Fen/Mgn ratios can be found up to 
redshifts of z ~ 6 (Freudling et al. 2003). However, 
the scatter of the Fen/Mgn ratios at z > 4 given 
in these studies is quite large. We presume that a 
significant fraction of the large range of these mea- 
sured Fen/Mgn ratios is caused by the restricted 
wavelength coverage of the analyzed quasar spec- 
tra. Although the mathematically best solution 



was determined for the multi-component fit, the 
continuum strength appears in some cases not well 
defined, yielding quite uncertain Fen flux mea- 
surements. Some scatter in the Feii/Mgn ratio 
given in recent studies may also be introduced by 
using emission lines that arc not completely cov- 
ered by the available spectra or measured at all, 
instead has been estimated from the Cm] A 1909 
emission blend (Freudling et al. 2003). 

In general, at lower redshifts (z < 2) the 
Fell/Mgn ratios derived from the SDSS quasar 
spectra (2.6 ± 0.3) are in good agreement with the 
ratios we measured (2.1 ± 0.3) from the composite 
spectra of z < 1 quasars with lower luminosities 
(see Fig. 4, luminosity sample, log\ La(1450A) < 
45 ergs -1 ). The sample of quasars studied by 
Wills et al. (1985) with z < 0.8 infer Fen/Mgll 
ratios which also are quite consistent with the 
ratios we derived. However, the intercalibrated 
Fell/Mgn ratios based on the SDSS quasars are 
significantly lower than the ratios we determined 
from composite spectra of nearly constant lumi- 
nosity for redshifts less than z ~ 2 (Fig. 5). A 
possible explanation may be given by the lower lu- 
minosity of these SDSS quasars compared to our 
composite quasar spectra in this redshift range 
with ; Off AL A (1450A) ~ 46.7 ergs -1 for the fol- 
lowing reasons. Although the MgnA2798 and Fell 
emission may both originate in the partially ion- 
ized zone of the line emitting region, both lines 
depend in different ways on the conditions of 
this region, e.g., radiative transfer, destruction by 
Balmcr continuum absorption, and conversion to 
optical Fen emission (e.g., Wills et al. 1985; Joly 
1987; Verner et al. 1999). Luminosity dependent 
variations in the SED of the AGN continuum, in 
particular of the ionizing continuum, are known to 
exist based on observational and theoretical stud- 
ies (e.g., Netzer et al. 1992; Laor et al. 1995; Wan- 
del 1999; Dietrich et al. 2002b). Hence, quasars 
of different luminosity may have a different con- 
tinuum SED which may result in slightly different 
conditions of the Fen and MgnA2798 emitting re- 
gion (z < 2). Taking the sensitive dependence of 
Fen and MgnA2798 on the conditions of the gas 
into account this may result in a wide range of pos- 
sible Fen/Mgn ratios as observed at low redshifts 
(z < 0.2). The large scatter of the Fe Ii/Mgll ratio 
also may emphasize that it is important to com- 
pare this line ratio for quasars of similar luminos- 
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ity. 

5. Discussion 

The Feli/Mgn emission line ratios which we 
find for our high-redshift quasar sample, as well 
as for the composite quasar spectra are consis- 
tent with results presented by Wills et al. (1985), 
Thompson et al. (1999), Murayama et al. (1999), 
Dietrich et al. (2002a), Iwamuro et al. (2002), 
and Freudling et al. (2003). In particular, the 
Feli/Mgn ratio at high redshifts shows similar 
strength compared to those measured in low red- 
shift quasars (z ~ 0.2) of comparable luminos- 
ity (WNW85). Assuming that the physical con- 
ditions in the BELR of quasars and that the spec- 
tral energy distribution of the ionizing continuum 
are independent of redshift the observed lack of 
evolution of the Fe Ii/Mgn emission ratio suggests 
that the observed relative abundance of iron, i.e. 
Fe/Mg, in the local universe is already achieved 
at redshifts z > 4. This result indicates that an 
episode of intense star formation began before the 
observed quasar activity. Hence, the nearly con- 
stant Fell/Mgll ratio over cosmic time has im- 
portant implications on the star formation his- 
tory, particularly at high redshifts since the age 
of the universe is less than ~ 1.3 Gyr (for redshifts 
z > 4.5). 

With this result, we can estimate the begin- 
ning of the first star formation in luminous quasar 
host galaxies, i.e., probably young massive ellipti- 
cal galaxies from the evolution time scale of type 
la SNe and the resulting delay in the Fe enrich- 
ment. Evolutionary models of the Fe enrichment 
in a galaxy following the initial starburst has been 
calculated e.g. by Yoshii et al. (1998). Although, 
as pointed out by Yoshii et al., evidence that 
the initial mass function may have been different 
for the first stars makes such model calculations 
quantitatively somewhat uncertain, it is clear from 
these computations that the relative Fe abundance 
is rather low initially and starts to grow steeply 
about <~ 1 Gyr after the beginning of the star for- 
mation, reaching a maximum (at about 3 Gyrs 
in the Yoshii et al. model), before declining to 
the local value (sec also Mattcucci & Padovani 
1993; Hamann & Ferland 1993). The same tem- 
poral evolution of the Fe/Mg ratio is found for the 
giant elliptical galaxy model (M4a) presented by 



Hamann & Ferland (1993). Their model predicts 
a strong increase of Fe/H at ~ 1 Gyr after the be- 
ginning of the star formation, with the most rapid 
rise of the Mg/Fc ratio during the 1 - 1.38 Gyr 
period. 

However, Matteucci (1994), Friaca & Terlevich 
(1998), and Matteucci & Recchi (2001) pointed 
out that the common assumption Td e i ay — 1 Gyr 
as a typical time scale for the enrichment delay 
of Fe, released by SNe type la SNe explosions, is 
only a rough estimate which based on the condi- 
tions of the solar neighborhood. The evolutionary 
models presented by Friaca & Terlevich (1998) and 
Matteucci & Recchi (2001) suggest that this delay 
depends crucially on environmental conditions of 
the star formation process and that it can be sig- 
nificantly shorter. They showed that for elliptical 
galaxies an intense early star formation episode 
with enhanced star formation rate and short star 
formation times scale results in a high SNIa rate 
which reaches a maximum already ~ 0.3 — 0.8 Gyr 
after the beginning of the star formation, in con- 
trast to <~ 4 to 5 Gyr as for the solar neighborhood. 
Various theoretical studies have explored possible 
scenarios that SNe type la can, under certain con- 
ditions, also be produced by massive stars with 
a much shorter evolutionary time scale than that 
assumed in the standard scenarios (e.g., Yoshii et 
al. 1998), some SNe type la progenitors could have 
much shorter life times (e.g., Iben & Tutukov 1984; 
Smecker-Hane & Wyse 1992). 

There is growing evidence that luminous quasars 
are residing in massive host galaxies (e.g., Kukula 
et al. 2001; Nolan et al. 2001; Dunlop et al. 2003), 
which at high redshift (z > 3) are probably young 
ellipticals or massive spheroidal systems in the 
process of forming (e.g., van Brcugcl et al. 1998; 
Papadopoulos et al. 2000; DeBreuck et al. 2002). 
This would support the assumption of significantly 
shorter evolutionary time scales for the major iron 
enrichment in massive spheroidal galaxies than 
that generally cited, r evo i ~ 1 Gyr. As noted 
above, the significant release of iron by type la 
SNe to the interstellar medium can be as short 
as r evo i < 0.3 — 0.8 Gyr. In addition, a more 
rapid evolution of progenitor stars to produce and 
release iron to the interstellar environment than 
previously assumed may be also possible. Re- 
cent models of galaxy formation, particularly in 
the early universe suggest that a rapid assem- 
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bly of (at least some) massive spheroidal systems 
was accompanied by intense star formation, and 
subsequent quasars activity (e.g., Kauffmann & 
Haehnelt 2000; Green et al. 2001; Romano et al. 
2002). The indication of an early star formation 
epoch beginning at Zf ~ 6 — 9 is also consistent 
with cosmic structure formation models (Gncdin 
& Ostriker 1997) and with recent estimates of the 
epoch of reionization. 

The observed lack of an evolving Fe n/Mgll ra- 
tio indicates that the star formation responsible 
for the chemical enrichment of the gas close to the 
quasar has started at least T evo i earlier. There- 
fore, the first major star formation epoch in high- 
redshift quasars (z > 4.5) started at Zf ~ 6 — 9. It 
is interesting to note that this estimate is compa- 
rable to the epoch of re-ionization of the universe 
(?Haiman & Loeb 1998; Becker at al. 2001; Fan et 
al. 2002). The indication of an early intense star 
formation epoch is further supported by recent 
results of the Wilkinson Microwave Anisotropy 
Probe (WMAP) (Bennett et al. 2003). Ccn (2003) 
suggests that the first stars were formed as early 
as z > 15 and that they re-ionize the universe for 
the first time, while a second re-ionization episode 
of the universe occurred at z ~ 6. 

6. Conclusion 

We have observed 6 luminous high-rcdshift 
quasars (> 4.4) in the near infrared wavelength 
range to measure the Fell/Mgll emission line ra- 
tio and thereby estimate the beginning of the first 
intense star formation. To explore the evolution of 
the Fen/Mgn ratio over a wide redshift range, i.e. 
cosmic time, we measured this ratio for compos- 
ite quasar spectra which cover a redshift range of 
z % 5 with nearly constant continuum lumi- 
nosity (redshift sample), as well as for composite 
quasar spectra which span nearly ~ 6 orders of 
magnitude in luminosity (luminosity sample). We 
find that the Fen/Mgil ratio shows no evolution 
for cosmic times corresponding 3.2 < z < 5.0 com- 
pared to measurements for quasars at z =~ 3.4 
(see e.g., Dietrich et al. 2002a). The Fell/Mgll 
ratio also shows no evolution from the local uni- 
verse to very high rcdshifts (z ~ 5), corresponding 
to an age of the universe of ~ 1 Gyr. For increas- 
ing continuum luminosities, however, we detect a 
weak trend for higher Fen/Mgn ratio at higher 



luminosities. It appears that at lower luminosi- 
ties (log\L x (U50) < lO^ergs- 1 ) Fell/Mgll~ 2 
while at higher luminosities (log ALa(1450) > 
lO^ergs" 1 ) Fell/Mgll~4-5. 

Assuming that luminous quasars reside in 
massive galaxies, probably elliptical galaxies or 
in forming massive spheroidal systems in the 
early universe (e.g., Kauffmann & Haehnelt 2000; 
Kukula et al. 2001; Nolan et al. 2001; Romano et 
al. 2002; Dunlop et al. 2003) the lack of evolution 
in the observed Fell/Mgll ratios across the red- 
shift range < z < 5 suggests that there is also no 
evolution in Fe/Mg in the quasar nuclear region. 
The actual value of Fe/Mg in quasar BELRs re- 
mains uncertain. Further modeling of the Fell 
and Mgn line emission is required. However, 
there is some evidence that Fe/Mg is above so- 
lar in quasars at low to moderate redshifts (Wills 
et al. 1985). If that result is confirmed, the lack of 
evolution in Fell/Mgll out to z ~ 5 reported here 
would indicate that the first major star forma- 
tion in the environment of high-redshift quasars 
started at Zf ~ 6 — 9, corresponding to an age of 
the universe of ~ 0.5 Gyr. It is interesting to note 
that this epoch is comparable to the proposed 
time of the re-ionization of the universe (Haiman 
& Loeb 1998; Cen 2003). 
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Fig. 1. — The restframe quasar spectra to- 
gether with results of the multi-component analy- 
sis are shown for BRI 0019-1522, BR 0103+0032, 
PSSJ 0248+1802, PC 1158+4635, Q 2050-358, 
and BRI 2237—0607. In the top panel the quasar 
spectrum is shown together with the power law 
continuum fit (dotted line), the scaled and broad- 
ened Fe-emission template, the scaled Balmer con- 
tinuum emission, and the Gaussian components to 
fit the MgnA2798 emission line profiles. The re- 
sulting fit is over plotted as a thick solid line. In 
the bottom panel the quasar spectrum is shown 
after subtraction of these components. 

Fig. 2.— The spectrum of SDSS 1204-0021, ob- 
served in the H-band spectral region, transformed 
to the restframe (top panel). In the middle panel 
the mean z ~ 4.5 quasar spectrum is shown, 
scaled to match the overall spectral shape of 
SDSS 1204-0021 and the resulting difference spec- 
trum is plotted in the bottom panel. 

Fig. 3. — The Fen/Mgn line ratio of the z ~ 4.5 
quasars in comparison with quasars at z ~ 3.4 
which have been observed previously (Dietrich et 
al. 2002a) are plotted as filled diamonds as a func- 
tion of redshift. The Fcii/Mgn ratio of composite 
spectra (filled boxes, redshift sample) and compa- 
rable luminosity like the individual high-redshift 
quasars are displayed to illustrate the lack of evo- 
lution of the Fell/Mgll line ratio. 

Fig. 4. — The same as Figure 3 but the Fe ii/Mgn 
ratios of the individual high-redshift quasars (filled 
diamonds) and composite spectra with increasing 
luminosity (filled boxes, luminosity sample) are 
plotted as a function of intrinsic continuum lumi- 
nosity, AL A (1450A). 

Fig. 5. — The same as Figure 3 but the Fe Il/Mgll 
ratios of the z — 3.4 and z — 4.5 high-redshift 
quasars (filled diamonds ♦) and redshift sample 
quasar spectra (filled boxes ■) of our study are 
compared to the results of Wills et al. (1985) □, 
Thompson et al. (1999) *, Iwamuro et al. (2002) 
(o, A), and Freudling et al. (2003) (small □). 
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Table 1 

The high redshift quasars sample 



quasar 


site 


date 


exposure time 








J 


H 

[sec] 


K 


BRI 0019-1522 


CTIO 


16-Scp-00 


2100 


2100 


2100 


BRI 0103+0032 


CTIO 


16,17-Sep-00 


8700 


8700 


8700 


PSS J0248+1802 


CTIO 


17-Sep-00 


7200 


7200 


7200 


PC 1158+4635 


Keck 


20-May-00 


2400 


2800 


1920 


SDSS 1204-0021 


VLT 


ll-Apr-01 
14-Jun-01 




1800 
1800 




Q 2050-359 


CTIO 


17-Scp-00 


5700 


5700 


5700 


BRI 2237-0607 


CTIO 


15,16,17-Sep-00 


15800 


15800 


15800 




12 3 4 5 42 44 46 48 

redshift z log XL X (1450 A) 
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Table 2 

MgiiA2798 and Feii emission line flux measurements of the z = 3.4, z = 4.5 quasar samples 

AND OF THE REDSHIFT AND LUMINOSITY SAMPLES. 



quasar 


z 


/o 5 AL A (1450A) 
[ergs" 1 ] 


F(Mgn) a 


UV Fcii/Mgn 


The z=3.4 Sample 


Q 0103-260 


3.375 


46.83 


73.6 + 3.5 


4.14 + 0.50 


Q 0105-2634 


3.488 


47.46 


267.4+ 11 


3.84 + 0.44 


Q 0256-0000 


3.377 


47.14 


200.7 + 7 


2.95 + 0.26 


Q 0302-0019 


3.286 


47.24 


341.5 + 10 


4.16 + 0.40 


Q 2050-359 


3.508 


47.43 


566.5 + 70 


4.10 + 0.80 


Q 2227-3928 


3.438 


47.07 


156.0 + 5 


3.90 ± 0.40 


Q 2348-4025 


3.310 


47.27 


155.6+10 


3.31 + 0.38 


The z=4.5 Sample 


BRI 0019-1522 


4.528 


46.92 


135.9+15 


4.15 + 0.66 


BR 0103+0032 


4.443 


46.99 


85.0 + 8 


4.58 + 0.67 


PSSJ 0248+1802 


4.443 


47.24 


134.3+13 


4.51 + 0.88 


PC 1158+4635 


4.733 


47.39 


394.8 + 40 


2.34 + 0.49 


BRI 2237-0607 


4.558 


47.27 


138.5+14 


5.42 + 0.95 


The Luminosity Sample 


comp.L401 


0.003 


42.50 


1003 ± 75 


2.43 + 0.34 


comp.L40h 


0.034 


43.51 


263 ± 10 


5.61 ± 1.23 


comp.L411 


0.051 


43.97 


273 ± 13 


1.68 + 0.40 


comp.L41h 


0.088 


44.46 


258 + 9 


2.23 + 0.37 


comp.L421 


0.267 


44.95 


105 + 3 


1.91 + 0.41 


comp.L42h 


0.721 


45.40 


38 + 2 


3.51 + 0.66 


comp.L431 


1.349 


45.94 


60 + 2 


3.53 + 0.86 


comp.L43h 


2.069 


46.45 


81 + 3 


4.43 + 0.94 


comp.L441 


2.780 


46.90 


131 + 4 


5.86+1.29 


comp.L44h 


2.909 


47.32 


312 + 21 


3.87 + 0.86 


The Rcdshift Sample 


comp.red025 


0.24 


46.29 


1278 ± 36 


4.54 + 0.43 


comp.red075 


0.80 


46.43 


249 + 6 


4.99 + 0.79 


eomp.redl25 


1.28 


46.56 


203 + 6 


5.45+1.08 


comp.redl75 


1.87 


46.66 


154 + 3 


4.02 + 0.53 


comp.red225 


2.14 


46.70 


122 + 7 


3.98 + 0.83 


comp.red275 


2.79 


46.73 






comp.red325 


3.17 


46.74 


72 + 5 


3.80 + 0.85 


comp.red375 


3.84 


46.82 






comp.red425 


4.26 


46.76 


64 + 5 


5.42 + 1.25 


comp.red475 


4.79 


46.65 


77 + 8 


2.23 + 0.39 



a [10 15 ergs 1 cm 2 ] 
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